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Abstract 

We study the possible indirect neutrino signal from dark matter annihilations inside the solar 
interior for relatively light dark matter masses in the 0(10) GeV range. Due to their excellent en- 
ergy reconstruction capabilities, we focus on the detection of this flux in liquid argon or magnetized 
iron calorimeter detectors, proposed for the next generation of far detectors of neutrino oscillation 
experiments and neutrino telescopes. The aim of the study is to probe the ability of these detec- 
tors to determine fundamental properties of the dark matter nature such as its mass or its relative 
annihilation branching fractions to different channels. We find that these detectors will be able to 
accurately measure the dark matter mass as long as the dark matter annihilations have a signifi- 
cant branching into the neutrino or at least the r channel. We have also discovered degeneracies 
between different dark matter masses and annihilation channels, where a hard r channel spectrum 
for a lower dark matter mass may mimic that of a softer quark channel spectrum for a larger dark 
matter mass. Finally, we discuss the sensitivity of the detectors to the different branching ratios 
and find that it is between one and two orders of magnitude better than the current bounds from 
those coming from analysis of Super-Kamiokande data. 
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I. INTRODUCTION 



The first exciting years of data from Large Hadron Collider (LHC) have provided us with 
the discovery of a Higgs boson that so far fits very accurately the properties expected from 
the Standard Model (SM) [l, 2]. Unfortunatelly, no significant deviation signaling physics 
beyond the SM has been found yet. Thus, the existence of dark matter (DM) remains one 
of our few confirmed signals of new physics. Sadly, all our present evidence for DM stems 
from its gravitational effects and thus we remain ignorant of its fundamental particle physics 
properties, such as its mass or couplings to other particles, that could help us understand 
the nature of DM. For this task, the LHC is still a very promising tool, as testified by the 
significant regions of the possible DM parameter space that it has already probed. In this 
work we will instead study a very complementary probe to the LHC searches which also 
has the potential to discriminate the DM mass and interactions: its indirect detection via 
annihilations in the solar interior giving rise to energetic neutrino fluxes js-Sj. 

For weakly interacting massive particles in the mass range between 0.1-1 TeV, such 



neutrinos could be detected by neutrino telescopes 
there have been several works 
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6|-[l0|, such as IceCube. In addition. 



16| on how well the same study would work for lighter DM 



with a mass of 0(10) GeV, and some bounds have been obtained for the Super-Kamiokande 



detector 



17H19| . These studies have focused on the prospects of running or planned detectors 



typically considered as far detectors of long baseline experiments. The detector technologies 
under consideration have been ranging from water Cerenkov detectors to liquid scintillator, 
magnetized iron (MIND) and liquid argon (LAr). 

A relatively common trait of the studies performed thus far is the assumption that, since 
the neutrino energy spectrum has a sharp cutoff at the DM mass, the DM mass will be 
very well determined by any experiment which is sensitive enough to detect the signal. 
However, two facts can hinder the DM mass reconstruction: (i) The neutrino energy will 
not be perfectly reconstructed by the experiment, (ii) If the neutrinos are not produced 
through direct annihilation into neutrinos, but rather through the decay products of other 
annihilation channels, then the neutrino spectrum will be softer and it will not always be 
easy to observe the cutoff energy. This is particularly true when there is a significant amount 
of background. The determination of the DM mass at the Deep Core experiment has been 
studied separately for the r and b channels 20|. The aim of the present study is to study 
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the simultaneous determination of the DM branching ratios and the DM mass at future 
neutrino detectors. 

The remainder of this work is organized as follows: In Sec. HIl we briefly review the 
physics of indirect neutrino signals from the Sun, starting from the gravitational trapping 
of DM inside the Sun, via the generation of neutrino events, to the neutrino propagation to 
and interaction in a detector on Earth. We continue by discussing the detector technologies 
used in the study in Sec. Illlt before we give the details of our numerical simulations in 
Sec. [IVl Section |V] is devoted to the presentation of the outcome of our simulations and the 
corresponding discussion. Finally, in Sec. IVTI we briefly summarize the results and present 
our conclusions. 



II. BACKGROUND PHYSICS 



Indirect neutrino signals from DM annihilations in the Sun and the Earth have been 



studied extensively in the literature (3|-|36| , in particular in connection to neutrino telescopes. 
The general idea is based on DM particles from the galactic DM halo scattering on the solar 
matter and losing enough energy to become gravitationally bound to the Sun. The DM 
then loses more and more energy in subsequent interactions before finally thermalizing in 
the center of the Sun, where they can annihilate with other DM particles into SM ones. Since 
neutrinos are the only SM particles that can escape the solar interior, these annihilations 
can be probed by searching for the resulting neutrinos at the end of the decay chain. 

The DM cross-section on SM matter can be effectively divided into a spin dependent (SD) 
and a spin independent (SI) cross-section. The latter of the two scales as A^, where A is the 
number of nucleons within the nucleus on which DM scatters, due to constructive interference 
among the scatters on individual nucleons. The SD cross-section, on the other hand, depends 
on the overall spin of the nucleus. Therefore, the bounds on the SI cross-sections from direct 
detection experiments tend to be much stronger than the corresponding bounds on SD cross- 
sections, since typically nuclei with a large number of nucleons are exploited. As the Sun 
consists mostly of elements with a small number of nucleons, the SD cross-section instead 
tends to dominate and it is mainly this type of cross-section that neutrino experiments will 
be able to probe. We therefore focus on the SD cross-section throughout this study. 

The typical capture rate of DM particles in the Sun by a species i can be described by 
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a = I dV ^B^a, fdu^ I dE\F{E)\\ (1) 



the relation [23l \3 



2mDMK J " 

where pdm is the local DM density, piir) is the density of the species i at a distance r from the 
Sun's center, m-QM is the dark matter mass, /ij the reduced mass of the DM-z system, f{u) 
the velocity distribution function for DM, which we will assume Maxwellian, and F{E) is a 
form factor for the scattering on i as a function of the nuclear recoil which is essentially 1 
for the dominant scatter in hydrogen at all relevant energies. For details on this procedure, 
see jiol . I37I 1 . For the purpose of our study, the main thing to notice is that the capture rate is 
directly proportional to the capture cross-section cTj. Taking only capture and annihilation 
into account, the number of DM particles in the Sun follows the differential equation 

^ = 5^ a - CaN' ^Cc- CaN\ (2) 

i 

where Ca is a constant describing the rate of annihilations between two DM particles. 
Since every annihilation involves two DM particles, the annihilation rate = CaA^^/2. 
In the steady state limit, where dN/dt = 0, we have detailed balance between capture and 
annihilation 

= ICc. (3) 

Thus, also the annihilation rate is directly proportional to the capture cross-section in this 
case. The the typical time-scale to reach the detailed balance is given by 

Tdb ~ J-^ (4) 

and the requirement that detailed balance is reached during the age of the Sun therefore 
translates to CcCa ^ (5 ■ 10^ years)"^. If Ca is fixed such that the DM abundance is 
provided as a thermal relic, then this puts a lower limit on Cc and thus also on a, the 
DM-nucleon cross-section. For the masses treated in this study, the detailed balance starts 
being in question for cross-sections cr ^ 5 ■ 10"'^ fb. As we will see, such small cross-sections 
cannot be probed by the detectors in question, regardless of the annihilation channel, and 
we therefore assume that detailed balance holds. 

Fixing Vq = 220 km/s, pdm = 0.3 g/cm^, and using the BP2000 solar model 
composition of the Sun, the annihilation rate F^ p into pp is thus given by 
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for the 



Ta,p = BR(DMDM ^ pp)F^ = BR(DMDM ^ pp)ag{mi^M) = aBrp^(mDM) (5) 



where we have introduced the parameters aBip = crBR(DMDM — )■ pp) and g{mDu) is a 
function of the DM mass which can be computed from the expression for the capture rate. 

In this study, we will be concerned with the annihilation rates into three different anni- 
hilation channels: 



DMDM — y qq, DMDM — y t't^, and DMDM — > vV. (6) 

Here, q denotes a generic charm or bottom quark channel, normalized to the annihilation into 
bottom quarks, while v denotes a generic annihilation directly into any flavor neutrino. As 
demonstrated in [l^, the neutrino spectra of the individual channels within these generic 
categories have negligible differences and the inclusion of separate degrees of freedom for 
each would only lead to almost perfect degeneracies among the components where only the 
following combinations can be independently constraint: 

(xBr^ = aBrfo + 0.4 crBrc and crBr^, = ctBt^^ + (jBt^^ + crBr,^^. (7) 

We now have a four-dimensional parameter space containing moM? O'Big, crBr,-, and 
crBrj,, all of which are related to the properties of DM. Depending on the parameters, the 
resulting neutrino spectrum may have very different characteristics, from the highly peaked 
mono-chromatic neutrino channel to the very soft quark channel, and for a varying range 
of different DM masses. We therefore pose the question of how sensitive future neutrino 
detectors will be to these parameters, since determining them would give us important 
insight to the properties and nature of DM. This is a subject that has already been partially 



studied in several contributions to the literature |12l . Il4j . |20| . ranging from huge neutrino 
telescopes to the smaller scale neutrino detectors used in terrestrial experiments. However, a 
common approach within the literature has been to fix some of these parameters, usually the 
DM mass, while exploring the effect of the others. In this study, we extend these previous 
analyses to the whole parameter space to explore possible parameter degeneracies as well as 
the capability of these detectors to actually measure the DM mass on their own (without 
fixing it). 

The most stringent upper bounds on the spin dependent cross-section of DM on protons 



in the mass region 10-25 GeV are given by the COUPP |39| and SIMPLE [40 



and are of (9(10) fb. The limits from an analysis of the Super-K bounds 



collaborations 



191 are around 



0.2 fb for annihilations directly into neutrinos, 1 fb for annihilations into r^r , and 10 fb 



for annihilations into bb (again fixing all other parameters). The values of the cross- sect ions 
that we will use as examples in our results section are all compatible with these limits. Note 
that, since a branching ratio can be at most one, if a signal is seen in the type of search 
discussed in this study, then the real capture cross-section must be as large or larger. At the 
same time, the cross-sections studied in 19| assume a 100 % branching into the respective 



channels and should therefore be treated as bounds on the very parameters that we are 
studying. Thus, crBrj, = 0.1 fb would not necessarily imply a large branching ratio into 
neutrinos, even though the quoted Super-K bound is 0.2 fb, since the actual capture rate 
may be as large as 10 fb. 



III. DETECTOR TECHNOLOGIES 



Large underground laboratories have been essential for the detection and study of very 
low luminosity processes such as neutrino physics or proton decay due to the necessary 
overburden to avoid additional backgrounds induced by atmospheric muons. The proposed 
next generation of big underground neutrino detectors aims to continue and build upon the 
success of present experiments entering the 100-1000 kt mass range. These detectors could 
play a key role in the future beam based oscillation experiments as well as serving as neutrino 
telescopes for astrophysics studies and increasing the sensitivity to proton decay. With the 
recent discovery of a large leptonic mixing angle ^13 by reactor neutrino experiments, this 
new generation of neutrino experiments could provide answers to fundamental questions 
about the neutrino mass hierarchy and leptonic CP-violation. With the increased angular 
and energy resolutions of these experiments, they will also provide a good opportunity for 
probing eventual neutrino fluxes from DM annihilations inside the Sun for DM masses in 
the (9(10) GeV range, which cannot be reached by neutrino telescopes due to a too large 
spacing between the optical modules. 

In this study, we focus on the liquid argon (LAr) time projection chamber and magnetized 
iron calorimiter (MIND) detector technologies. While water Cerenkov and liquid scintillator 
technologies could also be of interest, their capabilities for angular and energy resolution 
drop for energies above the quasi-elastic regime. For the detector technologies included 
in this study, we follow closely Ref. [l^, where the very same detector technologies were 
investigated. 
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a. Liquid Argon Time Projection Chamber: The LAr TPC is a modern version of 
the bubble chamber. By applying an electric field to very pure argon, ionizing particles 
passing through the detector can free electrons which subsequently drift to the instrumented 
surface of the detector where a three-dimensional read-out can be performed. This detector 
technology has been developed largely by R&D performed by the ICARUS collaboration 41l . 



42| and LBNE |43|, M] 



45|. 



It is under consideration as a detector for the future LAGUNA 
projects. LBNE is currently planned as a 10 kt detector on the surface in the initial stages 
of the project. However, the planed next stage of the project consists of an upgrade to 
34 kt underground detector at a later time. Similarly, the first stage of the LAGUNA LAr 
detector would correspond to 20 kt with the aim of a gradual upgrade up to 100 kt 
Here we will consider the physics potential of these upgrades assuming a 34/100 kt setup, 
keeping in mind that our results then scale in time, i.e., our assumed 10 years of running for 
the 34 kt detector would correspond to 34 years of running for the 10 kt detector (although 
this assumes that the detector is moved underground). Due to the difficulty in applying a 
magnetic field to the volume of liquid argon, we consider this detector type as being able 
to distinguish only the fiavors of neutrinos while having no sensitivity to whether or not an 
interaction was caused by a neutrino or anti-neutrino. 

b. Magnetized Iron Calorimeter: The MIND 46| technology consists of iron modules 
with instrumented layers in between, where the energy deposit of charged particles can be 
measured. Since electrons and positrons will not pass through the iron, this technology is 
only sensitive to muon-type neutrinos. On the other hand, a volume of iron can be eas- 
ily magnetized, providing sensitivity to distinguish neutrino from anti-neutrino events. A 
100 kt MIND detector is considered the baseline for implementation in a neutrino factory 
experiment. Currently, the India-based Neutrino Observatory (INO) detector is aiming 
for a MIND-type detector of 50 kt (alternatively 100 kt after upgrade) primarily in order to 
study atmospheric neutrinos but also as a possible far detector of a neutrino beam exper- 
iment. As for the LAr detectors, we will here study the capabilities of the fully upgraded 
100 kt version. The results would be equivalent to those for a 50 kt detector with twice the 
running time. 

For further details of our experimental setup, we refer the reader to [l^, in particular 
Tab. I and Fig. I, where the characteristics of the detectors are described in detail. 
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IV. NUMERICAL IMPLEMENTATION 



Our numerical implementation is an extension of that of [l^ in order to include the 
possibility of varying the DM mass within the simulations. We perform the following steps 
in order to accomplish this: 



1. Simulate the neutrino spectra. We start by simulatin 



three annihilation channels using t 



using the nusigma 



le Wimp Sim software 



49|, DarkSUSY 



50|, and Pythia 



; neutrino spectra for all 
which in turn is 



21 



48 



51| codes. While computing the 



spectra during the simulations for the specific masses involved would be preferable, this 
is too computationally expensive to be feasible and instead we simulate the spectra 
for masses between 5 and 30 GeV with a 0.5 GeV resolution and later interpolate 
between these once the results have been binned into 1 GeV bins. We chose 0.5 GeV 
in order to have a manageable amount of data while still having a resolution smaller 
than the bin size for better interpolation. For the direct annihilation into neutrinos, we 
assume an equal branching into the different neutrino channels in order to suppress 



oscillation effects on the monochromatic peak in this flux [21|, |33N35|. which could 
otherwise interfere with the results due to the solar oscillation length being of the 
same order of magnitude as the difference in the Sun-Earth distance between aphelion 
and perihelion for neutrinos with energies 0(20) GeV. The resulting fluxes would 
thus be strongly dependent to a change in the neutrino mass squared difference Am|^ 
within its current errors. For the other simulations, we used the neutrino oscillation 
parameters summarized in Tab. |T] in order to compute the fluxes. Since the quark 
and tau channels give rise to continuous fluxes rather than monochromatic ones, the 
energy resolutions of the detectors will give an averaged effect. It should be noted that 
the exact values of the neutrino oscillation parameters are not of crucial importance 



here, as oscillations main 



y act to equilibrate the fluxes among the neutrino flavors 



for continuous spectra (2l| . We therefore limit this study to a fixed set of oscillation 



parameters in the normal neutrino mass hierarchy. 



2. Simulate the detector response. 



input to the GLoBES software 



53 



'he fluxes from WimpSim were then used as 
54j in order to simulate the detector response. 



We bin the events in energy into 28 bins of size 1 GeV from 2 GeV to 30 GeV. The 
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Angles [°] 


Mass splittings [eV^] 


Phase 


012 


33.65 




7.54 • 10-5 


5 3.39 




38.41 




2.39 • 10^3 




613 


8.93 







TABLE I: The oscillation parameters used in our simulations. The mixing angles and mass square 
differences have been chosen to correspond to the best-fit values of 53]. Since the global errors are 



small for our purposes and the result is insensitive to the precise values of the parameters, we do 
not study variations of the parameters in this study. The same is true for changing from inverted 
to normal mass hierarchy. 

GLoBES flux normalization was cross-checked against a simple integration for exact 
energy resolution and 100 % efficiency in order to obtain the correct event rates. 
Finally, the detector resolution and efficiencies were switched on and the event rates 
were extracted for a capture cross-section of 1 fb and an exposure of 1 year. The 
detector response for different cross-sections and exposures were then obtained by 
simply scaling the response. 

3. Simulate the background. The main background to our signal processes is that 
of atmospheric neutrinos interacting in the detector. As atmospheric neutrinos come 



from all directions, this background can be reduced 
purpose, we use the atmospheric neutrino fluxes of 



Dy an angular cut. For this 
and oscillate them with the 



same oscillation parameters as those used for the signal simulation. For the same 
energy binning as in the previous step, we compute the atmospheric background in 
the detector. This is done by taking a sky average and then only including events 
within an angle from the Sun of the typical angular resolution of the detector. The 
assumptions for the angular resolution of each detector as a function of the energy 

n r 

are summarized in Fig. 1 of Ref. [1J|. We use the fluxes at the Frejus site from [55 1 
as representative when simulating the atmospheric background. As for the signal, the 
atmospheric background rate is computed for one year of running and later rescaled 
to the appropriate running time. 
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4. Implementation into a Markov Chain Monte Carlo. By interpolating the 
binned spectra, we can now obtain an approximation for the expected spectrum at 
an arbitrary DM mass. We can thus perform a Markov Chain Monte Carlo (MCMC) 
scan of the parameter space for any given parameter values within our restricted 
mass range. We perform simulations for the detectors outlined in the previous sec- 
tion and varying assumptions on the branching ratios and rriDu- The MCMC scan is 



implemented through the MonteCUBES software |56|], which employs the Metropolis- 
Hastings algorithm in order to sample the posterior distribution. This sampling of 
the posterior distribution is then used to construct credible regions through the use 
of the MonteCUBES Matlab GUI. For each MCMC scan, we have generated a total 
of four chains with 2 ■ 10^ samples each in order to have a reasonably well sampled 
posterior distribution. The convergence was good in all of the simulations presented 
in this study {R < 1.001 [s^j). Our MCMC simulations always assume a running time 
of 10 full years. 



V. RESULTS 



A. Detector rates 

In Fig. m we present the full binned spectrum in the energy range 2-30 GeV for the 
different annihilation channels assuming a 100 kt detector and a 1 fb capture cross-section. 
As can be seen from this figure, the neutrino annihilation channel offers the sharpest peak 
in the spectrum and is therefore naturally expected to give the best sensitivity to the DM 
mass. Furthermore, the peak of the neutrino channel coincides with the DM mass. The r 
channel spectrum is broad and hard enough to be visible above the atmospheric background 
rate and can also be exploited to provide some sensitivity to the DM mass. Finally, the 
softer and weaker quark channel spectrum will have a harder time providing a good enough 
signal-to-background ratio anywhere in the spectrum and the DM mass determination using 
the spectral cut-off will be very difficult unless capture cross-sections larger than 1 fb are 
invoked. Furthermore, in all cases, the largest expected event rates reach a few per year, 
indicating that a full run of (9(10) years will be necessary to obtain reasonable statistics. 
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FIG. I: The expected muon event rates in the different detector types for the different annihilation 
channels under study. The left (right) panels correspond to a MIND (LAr) detector of mass 100 kt. 
The assumed mass of the DM particles is 10 GeV (25 GeV) in the upper (lower) panels. The fluxes 
are produced assuming aHvx = 1 fb and full branching into the channel. Since LAr detectors 
do not have charge identification, we show the sum of fi~ and events for the plots using LAr 
technology, while the plots for the MIND detectors are for fi~ events only. 

B. Mass determination 



In Fig. ini we show the reconstructed 90 % credible regions for the mean outcomes of 
the experiments for a DM mass of 25 GeV and different annihilation channels. For each 
channel, we show examples for three different values of aBr^, except for the case of anni- 
hilation into quarks, where we show only two values in order to keep the figure reasonably 
clean. As already anticipated, we find that the hardest channels provide the best DM mass 
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FIG. II: The 90 % credible regions for the mean experimental outcomes for annihilations of moM = 
25 GeV DM into neutrinos (upper left panel), taus (upper right panel), and quarks (bottom panel). 
In all panels, the thick (thin) solid lines correspond to the 100 kt (34 kt) LAr detector, while the 
dotted lines correspond to the 100 kt MIND detector. For the upper two panels, different lines for 
the same experiment correspond to different values for aUr^/.,-: aBr,^ = 0.05 (blue), 0.1 (green), 
and 0.15 fb (red) for the left panel and crBrr = 0.1 (blue), 0.4 (green), and 0.7 fb (red) for the right. 
Due to the significant overlap in the lower panel, we show the results only for aBr^ = 3 (blue) and 
5 fb (red). 



determination. As can be seen from the two upper panels of the figure, the high energy 
resolution for the LAr detectors typically give good results for determining the DM mass 
for the neutrino and tau channels, providing error bars of around 1 GeV for the tau channel 
which improve by almost an order of magnitude for the neutrino channel. For the worse 
energy resolution assumed for the MIND detector the reconstructed error bars increase, with 
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1 GeV for the neutrino channel with the largest cross section and presenting significantly 
worse mass determination compared to the LAr detectors in all cases. The increased mass 
sensitivity for increased aBi^ is also naturally expected, as a higher number of events can 
be accumulated. In the opposite end, when aBr^. is close to the sensitivity limit of the 
experiment, very little sensitivity to the mass is achieved, given the low statistics and poor 
signal-to-background ratio. 

The lower panel of Fig. |Tll corresponds to the case in which DM is annihilating purely 
into quarks. As expected, the softer spectrum significantly hinders the reconstruction of the 
DM mass. Moreover, the plots also show clear signs of degeneracy for the crBig channel at 
lower DM masses than the 25 GeV assumed, which further decreases the ability to accurately 
determine the mass. In order to understand the origin of this degeneracy, in Fig. Illll we show 
the credible regions for several of the projections for the case of crBr^ = 3 fb at the 34 kton 
LAr detector. While this is enough in order to detect that there is a signal, it is not clear what 
the signal consists of. In particular, a strong degeneracy appears around mDM — 14 GeV. 
The existence of this degeneracy was confirmed by using the built-in degeneracy finder 



see 



58[) of MonteCUBES for all our experimental setups and for those masses where the 



degenerate solution could also be found within the region of study (10-25 GeV). For the 
case presented in Fig. Illlt the best-fit point in the degenerate likelihood maximum is given 
by: moM = 14.2 GeV, cxBr, = 1.3 fb, (xBr^ = 0.15 fb, and crBr^ = 9.0 ■ 10"^ fb. In order 
to visualize the origin of the degeneracy, we show in Fig. IIVI the fiux due to aBiq = 3 fb 
for mDM = 25 GeV along with the different contributions to the flux from the degenerate 
solution and their sum. From this flgure it is clear that, in the energy range where the signals 
are well above the background level and at the same time have an event rate which could 
be accessible within a 10 year time frame, the aBr^ and crBrT- components of the degenerate 
solution add up to form a spectrum which is indeed very similar to that of the original input. 
For the degenerate solution, the part of the spectrum in the aBiq channel that is lost by 
going to a lower DM mass is compensated by the harder spectrum of the crBr^ channel. 
This corresponds very well to what we would expect from Fig. IIIIl where we can see that a 
signiflcant crBr^ is always needed for the degenerate solution. At the same time, we can see 
that, as we decrease the mass towards that of the degenerate solution, aBig decreases while 
crBrT- increases correspondingly. It should also be noted that for the lightest DM masses 
in the degenerate solution, aBi^, increases signiflcantly. This is due to the r spectrum not 
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being hard enough at the low DM mass to reproduce the original signal without involving 
also the v channel in the observable bins with the highest energy. 

A similar degeneracy can be observed between a pure cxBr,- signal at lower masses and a 
combined degenerate signal dominated by crBrg at higher masses. However, this degeneracy 
is not as strong as that for the pure aBr^ signal, as it is slightly more difficult to combine 
the different signals to reproduce the harder crBrT--like spectrum. 



C. Sensitivity to different branching ratios 

In Fig. |V] we show the sensitivity of the different detectors to the various annihilation 
channels as a function of the DM mass. In order to derive these bounds, we compute the 
90 % credible upper bound on the different branching ratios in absence of signal. In these 
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FIG. IV: The event rate from annihilations of DM with moM = 25 GeV and aBvg = 3 fb in a 34 kt 
LAr detector is shown (red dashed hne) together with the event rates from the degenerate solution 
m = 14 GeV, aBvg = 1.3 fb (red solid line), ciBr^ = 0.15 fb (blue solid line), and crBr^ = 9.0-10^'' fb 
(purple solid line), and their sum (green thick line). Also shown is the atmospheric background 
rate (black solid line), while the horizontal gray line corresponds to an average of one event per 
bin over a 10 year period. 



a:25GeV,aBr = 3 fb 



b: 14GeV,aBr =1.3fb 




simulations the DM mass is kept fixed to the corresponding value in the horizontal axis in 
order to obtain the sensitivity for that particular mass and to be comparable with bounds in 
other studies. All other parameters have been marginalized over. While we have seen that 
the determination of the DM mass can be a challenging task in most of the cases presented 
above, it may still be of interest to pose the question of how well the different branching ratios 
can be constraint if the DM mass is known from a different source, such as direct or collider 
DM searches. As expected, we see that for all cases the sensitivity becomes significantly 
worse at low DM masses, a result both of the signal appearing in a smaller energy range 
and the larger atmospheric background at low energies. For all channels, we can observe 
that the sensitivity of the 100 kt LAr detector would be roughly half an order of magnitude 
(a factor of ~ 3) better than that of the 100 kt MIND, with the largest difference in the 
neutrino channel, where the better energy resolution of the LAr detector at higher energies 
plays a decisive role. As expected, the sensitivity to the branching ratios becomes better the 
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FIG. V: The 90 % credible upper bounds for the aJirx parameters for the (upper left), r (upper 
right), and q (lower) channels under the assumption of a mean experimental outcome in absence 
of signal as a function of the DM mass mDMj which has been assumed to be known. 



harder the spectrum for the channel is, with the sensitivity to crBr,^ ranging between 10~^ 
to 10"^ fb, the sensitivity to aBij. between 10~^ and 10~^ fb, and the sensitivity to crBig 
between 0.3 and 5 fb. 



VI. SUMMARY AND CONCLUSIONS 



We have examined the neutrino signal from DM annihilations inside the interior of the Sun 
in different large projected neutrino detectors using the liquid argon (LAr) and magnetized 
iron (MIND) technologies. Such detectors are being discussed for the next generation of 
far detectors for neutrino beam oscillation experiments as well as for neutrino telescopes 
and, as we have shown, would also allow to probe the nature of DM. The energy reach 
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of such detectors limits our study to relatively low DM masses 0(10) GeV, which is quite 
complementary to similar studies in neutrino telescopes such as IceCube that operate at 
higher energies (9(100 — 1000) GeV, but have limited sensitivity at lower energies due to a 
large distance between the optical modules. Our results are based on the detectors taking 
data for 10 full years. 

In this study, we have extended the work in the existing literature by allowing the DM 
mass to take varying values when performing our fits. Studies similar to ours have been 
presented before in the literature, but typically fixing the DM mass at particular values. 
Thus, we have been able to study how well the DM mass, along with the branching ratios 
of the annihilation, can be determined from the shape of the neutrino spectrum. 

We have shown that the best annihilation channel for a precise determination of the 
DM mass is the direct annihilation of DM into neutrinos, presenting us with an essentially 
monochromatic spectrum. The energy resolution of the detector is then of crucial importance 
for how well the DM mass can be determined, with error bars of (9(0.1) GeV for the LAr 
type detectors and a few GeV for MIND (at a DM mass of 25 GeV). For the relatively hard 
spectrum of DM annihilations into r^r+ pairs, we still have a fair chance of constraining 
the DM mass, at least for the cases where the branching ratio is sufficiently large. The LAr 
detectors now show a mass resolution of a few GeV at best (again for a DM mass of 25 GeV), 
while the MIND is again significantly worse. The worst sensitivity to the DM mass is found 
for the annihilation into quarks, which produces a very soft spectrum of neutrinos. Not only 
would this channel produce less neutrinos in total per annihilation, but it would also be 
challenged by a higher atmospheric background at low energies. In addition, the spectrum 
is so soft that the maximum energy bins for which events could be expected over a 10 year 
period (a rate larger than 0.1 events/year) do not correspond to the kinetic cutoff at the 
DM mass. The result is that the resolution for the mass suffers even for the LAr detectors. 

In addition to these effects, we have also shown the existence of a degeneracy in the 
determination of the DM mass and the branching ratios. This degeneracy originates in the 
fact that the harder spectrum for annihilations into r^r+ pairs for a lower DM mass can 
be quite similar to that of annihilations into quarks for a higher DM mass. While this 
does not pose a real problem for the actual indirect detection of the neutrino signal from 
DM annihilations in the Sun (there is a clear signal over the atmospheric background), it 
would do so for the measurement of the DM mass and branching ratios and therefore for the 
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determination of the nature of DM. It should be noted that if the DM mass (or branching 
ratios) would be determined from a different source, this degeneracy would be lifted. 

In conclusion, the large future neutrino detectors intended as far detectors of neutrino 
beam experiments also present an exciting opportunity to learn more about the properties 
of DM. In particular, the DM mass can be probed at these detectors, but its accurate 
determination strongly depends on the actual branching ratios of DM annihilations and, 
when all parameters are simultaneously taken into account, strong degeneracies may appear. 
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